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ABSTRACT 


A radiant heat transfer computer program has been developed to cal- 
culate radiation from inhomogeneous gases prevalent in rocket exhaust 
plumes from clustered engines. The infrared spectral absorption charac- 
teristics of the radiating species considered in this computer program - 
water vapor, carbon dioxide, carbon monoxide and carbon particles - have 
been determined. Band model parameters have been used to represent the 
infrared spectral absorption coefficients over 25 cm"^ increments. A 
modified Cur tis-Godson approximation, used in the inhomogeneous heat 
transfer calculation, has been shown to give satisfactory results over 
the temperature and pressure range of interest in Saturn exhaust plumes. 
Results are shown for the Saturn-type engines for specific flow field 
assumptions. Some comparison with experimental spectroscopic data will 
also be presented. The effects of wavelength increment, field of view, 
and distance increment along the line of sight on the heat transfer, as 
well as the computer techniques for minimum computer time in calculating 
radiation from a three-dimensional flow field, are discussed. 
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Definition 
spectral absorptance 

average absorptance over a wave number interval 
Planck function 

mole fraction of absorbing gas 
spectral line spacing 
band width 
optical path 

optical path for coll is ion-broadened 1 ines 
optical path for Doppler -broadened lines 
Ladenberg-Reiche function 
spectral absorption coefficient 
carbon spectral absorption coefficient 
path length 

hypothetical path in the Cur tis-Godson approximation 
temperature corrected pressure 
total pressure 

pressure of absorbing gas in the (homogeneous) h^^ zone 

pressure of foreign gas broadener in the (inhomogeneous) 
h^^ zone 

intensity distribution function of a band radiant flux 
per unit area 

radiant flux per unit area 
spectral line intensity 
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DEFINITION OF SYMBOLS (Continued) 

Definition 

average line intensity over a finite wave number 
increment 

hypothetical intensity from an inhomogeneous gas for 
spectral lines 

spectral line intensity per unit pressure 
temperature in degrees Kelvin 
reference temperature in degrees Kelvin 
mass of gas per unit area 
equivalent width 

argument of Ladenberg-Reiche function 

argument of Ladenberg function in the Curtis -Gods on 
approximation 

carbon mass fraction 

half-width of spectral line 

half-width of absorbing gas for homogeneous path 

half -width of foreign gas for homogeneous path 
band model collision-broadened half-width 
band model of Doppler-broadened half-width 
wave length 
wave number 

angular spherical coordinates 

ratio of half-widths due to self-broadening 

ratio of foreign gas half-width to water half-width 
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DEFINITION OF SYMBOLS (Continued) 


Symbol 

(7/d) 


H2O 


t(v) 


Definition 

mean half-width to line spacing for water 
spectral transmittance 
average transmittance 

inhomogeneous gas transmittance in the Curtis-Godson 
approximation 
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SUMMARY 


A radiant heat transfer computer program has been developed to cal- 
culate radiation from inhomogeneous gases prevalent in rocket exhaust 
plumes from clustered engines. The infrared spectral absorption charac- 
teristics of the radiating species considered in this computer program - 
water vapor, carbon dioxide, carbon monoxide and carbon particles - have 
been determined. Band model parameters have been used to represent the 
infrared spectral absorption coefficients over 25 cm“^ increments. A 
modified Cur tis-Godson approximation, used in the inhomogeneous heat 
transfer calculation, has been shown to give satisfactory results over 
the temperature and pressure range of interest in Saturn exhaust plumes. 
Results are shown for the Saturn-type engines for specific flow field 
assumptions. Some comparison with experimental spectroscopic data will 
also be presented. The effects of wavelength increment, field of view, 
and distance increment along the line of sight on the heat transfer, as 
well as the computer techniques for minimum computer time in calculating 
radiation from a three-dimensional flow field, are discussed. 


I. INTRODUCTION 


Heat transfer to the base of a launch vehicle can occur in two 
ways: by convection and radiation. Both of these modes are important 
especially on large vehicles which burn RP-1 (kerosene) and oxygen and 
on vehicles which burn hydrogen and oxygen. There is currently no com- 
pletely satisfactory technique of calculating base heating for a given 
engine array at various altitudes and vehicle velocities. Therefore, 
both a better means of measuring heat transfer to the base region of 
vehicles actually under test and a more rigorous and reliable means of 
predicting heat transfer to vehicles before testing are required. 

Convection to the base region depends upon the properties of the 
gas flow close to the surface being heated. For this type of heating, 
simple models of the flow, together with boundary layer correlations 
for the heat transfer coefficient, have been used in the past to arrive 
at local heating rates on the heat shield. Radiation, on the other 
hand, is sensitive to the entire flow field. Predictions have been 
made either by representing the exhaust plume as a surface with a given 
emissivity, temperature, and configuration, or by assuming a simple 



geometrical configuration for the hot gas, which is presumed to be uni- 
form, and integrating the equation of transfer with a simplified 
absorption coefficient. 

The exhaust plume from a multiple-engine configuration like the 
Saturn V is an extremely complex gas-dynamical structure. The exhaust 
jets from the clustered engines interact with the atmosphere, and with 
each other, through an intricate series of shock waves and turbulent- 
mixing layers (figure 1). At the higher altitudes, on a multi-engined 
vehicle, the high- temperature, high-pressure collision zones between the 
engines can be a strong source of radiation (figures 2 and 3). Wherever 
the exhaust material that flows back into the base region mixes with 
the atmosphere, burning will occur; at low altitude, this burning is 
substantial. With such a complicated flow field controlling the heat 
transfer to a base region, it is not surprising that simple models of 
the heat transfer must be tied strongly to empirical data. 

Helpful discussions with Mr, Werner K. Dahm on the radiative heat 
transfer problem and in the program definition are greatly appreciated. 
The cooperation of Mr, Robert Yossa, Mr. John Reardon and Mr. Marcus 
Dash is gratefully acknowledged, and assistance of Mr, Addison McGarrety 
in performing the computer programming is greatly appreciated. 


II, OUTLINE OF THE RADIATION PROGRAM 


In an effort to learn more about the base heating problem and to 
provide both mathematical techniques and values of parameters required 
for solution of the problem, a number of studies have been made both 
under contract and in-house at the Marshall Space Flight Center. These 
s tudies include the following; 

(1) Determination of accurate flow fields, i.e,, the tempera- 
ture, pressure and concentration field of the exhaust plumes for single 
and multiple rocket engines. 

(2) Accurate absorption coefficient data for the radiating 
species. This involves the determination of the band model parameters 
for the radiating species for the temperatures, pressures and concentra- 
tion typical of Saturn V exhaust plumes. Water vapor, carbon dioxide, 
and solid carbon particles (soot) are the principal radiators in the 
first-stage plumes. Water vapor is the only radiator of importance in 
the plumes of the hydrogen-fueled stages. Extensive experimental and 
theoretical studies of the emission of CO^ under various conditions have 
been made, and the properties of this molecule are well understood. 

Much less is known about H 2 O. 
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(3) An exhaust plume radiation calculation program for 
inhomogeneous gases . 

(4) Comparisons of calculated values with full scale radia- 
tion data. 


III. MSFC PLUME FLOW FIELD PROGRAM 

1 


For first stage vehicles, the heating environment is determined 
from total and spectral radiation measurements in single engine exhaust 
plumes, and an empirical model developed, using a total emissivity in 
conjunction with a form factor program to calculate the radiation to the 
base. Because of the contribution of the afterburning layers, this cal- 
culated heat flux is at a maximum at sea level. This maximum value is 
used in the design of the heat shield. The plume radiation decreases as 
altitude increases, and the plume spreads, thus becoming cooler and less 
dense. 

The upper stages of the Saturn V are powered by LH2*L0X engines, 
the Rocketdyne 200,000 lb thrust J-2. The S-II stage has a cluster of 
five J-2*s, while the S-IVB has only one J-2. The radiation from the 
first stage engines is principally carbon, and is therefore continuum 
radiation; but the principal radiator in the upper stages is water vapor, 
with a band spectrum. The first stage design information obtained from 
experimental data on single engines is a relatively simple problem, but 
the design information on the upper stages must be entirely theoretical 
since scaling laws for model engine data are not yet known. 

The chemistry of the LH 2 -LOX plume can be described if the appro- 
priate thermodynamic properties and velocity can be solved simultaneously. 
The plume flow properties of the single engine on the S-IVB stage can be 
predicted with an axisymmetric method-of-character is tics program. To 
obtain flow properties in the intersecting plume region of the S-II stage, 
a three-dimensional program is needed. The gas plume in the intersection 
region radiates heat much more strongly than the undisturbed plume since 
the gas in the intersection region has crossed a shock wave, causing its 
temperature and pressure to rise. There are presently two efforts being 
funded which will result in programs to predict the flow in the plumes 
and intersection region of the S-II. One study uses a three-dimensional 
method-of-characteris tics scheme. This scheme has the disadvantage of 
not being able to handle the imbedded subsonic regions which sometimes 
occur in the collision zones [2]. The other study uses a finite dif- 
ference technique [4], which should be able to handle both subsonic 
and supersonic flow. All of the three-dimensional flow field work now 
being done is restricted to inviscid, non-reacting, constant gamma flows. 
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Equilibrium flow capability may be added later. The free shear layer is 
being neglected in the computation of radiation from the J-2 plumes. 

The plume flow field work now being done by MSEC is summarized below. 

The items marked with asterisks are actually being used to furnish input 
to the radiation program; the other items should result in a better under- 
standing of plume flow phenomena and will be used in other studies, 

A. Single Nozzle 

1. Inviscid Core 

*a. Axisymmetric method-of-character is tics with either 
ideal gas or equilibrium. 

b. One-dimensional gas-dynamic flow with equilibrium or 
finite rate chemistry for H 2 -O 2 systems. 

2. Nozzle boundary layer with mass injection and either ideal 
gas or equilibrium chemistry. 

3. Turbulent free-shear layer with equilibrium chemistry for 
almost any gas, but finite rate chemistry for H 2 -O 2 systems. 

4. Multiple shock program. • 

5. Viscous far wake. 

B. Clustered Nozzles * 


*1. Inviscid three-dimensional method-of-character is tics with 
ideal gas [2] , 

*2. Inviscid numerical finite difference scheme with ideal gas. 


IV. BAND MODELS 


If, in a specific wavelength region, the physical state of an 
absorbing or emitting gas is known, as well as the locations, intensities, 
and shapes of the lines, it is possible to calculate the radiation emitted 
by a specific sample of the gas. However, for the gases of interest in 
Saturn- type exhaust plumes, the absorption coefficient varies extremely 
rapidly with wavelength. Consequently, an exact calculation of the 
radiant heat transfer would be a formidable task even with modern 
computers. By using band models, it is possible to replace this detailed 
calculation over frequency by an average over selected frequency intervals. 
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Another objective of the use of a band model concept is to permit the 
interchange, in the heat transfer equation, of the order of integration 
over spatial and frequency variables. Band models may be said to have 
their origin in the theory advanced by Elsasser in 1938. For band 
models, in general, one begins by considering the equivalent width of 
a single line, defined by 





( 1 ) 


where A is the absorptance. Assuming the Lorentz line shape, we have 
for the absorptance of a mass of gas, U, cm"^ (= Pi) 


A = 1 - exp[-K(v)U] , 


( 2 ) 


where 


K(v) = ^ 

(v - Vq)^ + 7 ^ 


(3) 


S is the spectral line intensity defined by 


s 



dv. 


o 


( 4 ) 


and 7 is the width of the line at ha If -maximum intensity. Integrating 
over all frequencies gives 


W = 2nyX e'^[I^(X) + Ii(X)], 


(5) 


where Io(X) and Ii(X) are pure imaginary Bessel functions and X = SU/2jt7, 
The importance of the Ladenburg-Reiche solution - i,e., W = 2Tiyf(X) where 
f(X) = X6 "^[Iq(X) + Ii(X)] - is that the absorptance determined from the 
statistical band model can be expressed in terms of the characteristics 
of a single line. Also, any curve of growth (that is, a graph of W versus 
U) must have the shape of this function. By appropriate application of 
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the experimental data and the theoretical curve of growth, the param- 
eters S/d and 7/d can be determined for a given pressure and path 
length. 


Two useful models are the Elsasser and the random models. The 
latter was proposed by Mayer [ 5 ] and by Goody [6]. Sometimes the random 
Elsasser model is used, but this has been shown to reduce, ultimately, 
to the random statistical model. 

A. The Elsasser Model 


The Elsasser model assumes an infinite array of equally spaced, 
equally intense lines, all of which have the Lorentz shape. Since the 
array is infinite in extent, any line center can be chosen as the 
center of the complete array. Thus, at a frequency v, the absorption 
from a single line is 


Kvi , KV2 — — , ...,Kv , 

(v - Id)^ + 7^ (v - 2d)^ + 7^ ^ (v - nd)^ + 7^ 


and the absorption at frequency v from all the lines is then 


K 


V 


n=cxD 



n=-oo 


(s/n) 7 ^ 

(v - nd)^ + 7^ 


( 6 ) 


If we write W = Ad , since the absorptance over any portion of 
the band is the same as the absorption over a line, the formula 


A = 1 - exp(-K^U), 


when integrated over dv, yields 




JL i 

r 

exp 

su 

s inh 

e 1 

2jt J 

d 

cosh p - 

cos 2 Z 


- 7 t 


( 7 ) 
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where Z = Jtv/d, and p = 2ziyld. Two limiting cases of this integral are, 
in the weak-line approximation, 


A(v) = 1 - exp(-px), 


and, in the strong-line approximation, 
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A(v) = f(it 7 SU/d^), 


( 8 ) 


where the right-hand side takes the form of an error function. 
B. Statistical Model 


The statistical model for a band assumes that the positions of 
the lines occur randomly and that the intensities of the lines in the 
band can be represented by some probability distribution function. The 
spectral transmittance, t, is written in the form of Beer*s law, and an 
average transmittance is written as 


J • J' J p(vi, V 2 , . . . , v^) dvi .. 


00 00 
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d V. 
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n P(Si)e dS. 


- _ Avi Avg Avj 


00 00 


p(Vi, . . . ,v. ) dVi,dV2 . 


. dv. 


“ P(S.) dS. 

jt 1 1 


Avi Av^ 
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The argument of the exponential function contains the integrated intens- 
ity of a single line whose center frequency is situated at v^. In terms 
of the absorptance, this expression reduces to 


But, since 



00 



W(S) P(S) dS 


o 


(9) 
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is the average value of the equivalent width for a single line, the 
absorptance is 


A(v) = 1 - exp[-W/d], 


( 10 ) 


where W/d as a function of pressure and path length can be determined 
from a curve of growth at a particular temperature and over the particular 
wave number interval. Curves of growth, W/ d versus U, constructed from 
various probability intensity distribution functions P(S) are shown in 
figure 4, and the corresponding absorptances are shown in table II. From 
figure 4, the maximum difference of W/d for the curves of growth using a 
delta, exponential and a 1/S probability intensity distribution function 
is approximately 25 percent. 


V. THE MODIFIED CURTIS -GODSON APPROXIMATION 


To calculate the radiant heating from rocket exhaust plumes, we 
must know the spectral transmittances of the exhaust plume radiating 
gases. Since exhaust plumes are strongly inhomogeneous, we must deter- 
mine the spectral transmittances of inhomogeneous gases in order to 
develop a calculation method. Under MSFC contract, a general method 
has now been developed to calculate spectral transmittances of inhomo- 
geneous gases from the properties of homogeneous gases. Thus, any 
spectral transmittance for a particular inhomogeneous gas path can be 
calculated by properly combining known data on gases at constant tem- 
perature, pressure, and concentration [7]. 

The method developed is based principally on two special spectro- 
scopic concepts: the molecular band model and the Cur tis- Godson approxi- 

mation. The band model yields an explicit closed formula for the 
molecular radiation within each selected spectral region of interest, 
which uses as input data the averaged line strength, the averaged line 
spacing, and the averaged line half-width. An average of 25 wave numbers 
is considered sufficient. The use of a band model is critically impor- 
tant for practical calculations of gas radiation. The Curtis-Godson 
approximation is a method of combining the parameters that appear in the 
band model formulae in such a way that the parameters needed for an 
inhomogeneous gas calculation are obtained solely from homogeneous gas 
data. This is a critical factor because the necessary parameters need 
be calculated or measured in the laboratory only for uniform gas samples. 
Without the Curtis-Godson or some equally good approximation, we would 
have to treat each inhomogeneous gas path as a special case. The band 
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model used was the random model with equal line intensities and equal 
line widths. In terms of the spectral transmittance, t, at wave 
number, v, 

lni = 2nJf(X), (11) 

where 


^ cs/d:)_ u 

2:t(7/d) - 


( 12 ) 


Since the Curtis-Godson method is based on the premise that one can 
substitute a hypothetical homogeneous path for an inhomogeneous path, 
this hypothetical path will have the same spectral transmittance as an 
inhomogeneous path if the parameters ( 7 /d)* and X* as defined below are 
used in the band model formulae for the transmittance of a homogeneous 
gas . 


If an inhomogeneous gas has a transmittance of t* , then some 
hypothetical homogeneous sample has the same transmittance if it has 
certain values of ( 7 /d)* and X*. The Curtis-Godson approximation gives 
us 


and 





-In T, 
1 


f(Xi) 




i 


2ir 


I 


X. 


-In 

1 


1 Lf(Xi> J 


(13) 


(14) 


which define the band model parameters of this particular homogeneous 
sample. The definition is given in terms of the band model parameters 
of the individual zones comprising the inhomogeneous path being studied; 
these zones are labeled with the subscript h. 
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-In t' 



(17) 


Generally, transmittances calculated with equation (17) are not very 
sensitive to errors in Therefore, with moderately accurate values 

of the accuracy with which t* can be calculated is limited by the 

accuracy of the transmittance measurements and the suitability of the 
theory. This is particularly true in the high and low X regions where 
equation (17) reduces to a form in which x’ is independent of Xi. Thus, 
if all the X values are low (less than 0.2), f(X) ~ X, and equation (17) 
reduces to 
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-In t' 


(18) 



(-In 



A 


If all the X^'s are high (greater than 2.0), and f(X) (2X/it)^^^, 
equation (17) reduces to 


“In T* 



(19) 


Equation (17), the general expression, can be used for any X values. 
Equations (18) and (19) are good in the low and high X regions, respec- 
tively. Experiments were carried out to test the Curtis-Godson theory 
for temperatures and pressures typical of rocket exhaust plumes. Band 
model parameters were measured for isothermal gases at various tempera- 
tures, the results were used to predict transmittances of known inhomo- 
geneous paths, and the inhomogeneous path transmittances were compared 
to the theoretical prediction. Using a furnace-gas cell arrangement, 
numerous transmittances of water vapor inhomogeneous paths were measured 
and compared to transmittances calculated by the Curtis-Godson combina- 
tion method in conjunction with the statistical band model. The results 
are shown in Tables III, IV, V, VI, and VII. 

VI. DETERMINATION OF THE BAND MODEL PARAMETERS 


Progress in the determination of the band model parameters and 
absorption characteristics for the radiating species in the Saturn-type 
exhaust plumes, H 2 O, COp, CO and carbon particles is described as 
follows : 


1. H 2 O 

From the band model formulae, one can see that if values of 
S/d(T,v) and 7 /d (T, v, g^g) can be determined, then the absorp- 
tance can be calculated for any P^. For thin-gas spectra where 
% (SU/d)/(47/d) « 1, S/d can be determined directly. A great many water 

vapor spectra over the wavelength range from 1 to 22 [j,, which have been 
published in recent years, cover the temperature range from 300°K to 
2700°K. The optical depths range between 0,2 and 100 cm-atm, and total 
^ pressures between 50 mm Hg and 10 atm. If sufficient thin-gas spectra 

at all temperatures were available, a set of S/d (T,v) values could be 
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derived directly. To obtain S/d (T,v) values from measured spectra in 
which the gas was not thin, the functional relationship between emissiv- 
ity and optical depth must be known. The curve of growth given by the 
statistical band model has been used to determine the mean spectral 
absorption coefficients from thin and non- thin spectra by adjusting the 
fine structure term in such a way that the integral of the absorption 
coefficients over a given vibration rotation band results in the known 
value of the band intensity. Their values have_been tabulated in refer- 
ence 3. In general, the mean deviation of the S/d (T,v) values is within 
±20 percent (figures 5 and 6 ). 

The average line spacing values, d(T), obtained from the 
value of 7 /d (T,v) using the expression 


2 

d 


(v) 


d(v) 



C(n/t^7t’ + a' ) + (1 - C) 
a 



where 


and 


a* 

a 


7° 


a 


b 
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( 20 ) 


and where C is the mole fraction of water vapor. The value of 7 ^ ^ is 
0.5 cm”^ atm"^ and a* = 0 . 1 , 2 . 

These results are shown in figure 7 along with the results 
of several other investigators. The straight line, a least squares fit, 
is given by 


d(T) = exp[-.00106T + 1.21], 


( 21 ) 




♦ 


tr 
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These expressions for S/d and 7 /d were compared with the data of Simmons 
Arnold and Smith [18]. (The results are shown in figure 8 .) The above 
mentioned values of S/d (T,v) and 7 /d (T.Pforeign gas) being 

used in the radiative heat transfer computer program. The values of 
7 /d for H 2 O now in use do not take into account the wave number depend- 
ence. Since the 7 /d values are determined from the asymptotic regions 
of the curves of growth [15], long path length measurements are required 
Consequently, a long burner program was initiated to determine the neces 
sary curves of growth. The design and fabrication of this burner have 
been completed under MSFC contract [3]. Measurements now in progress 
should be completed by November 1966. The long burner shown in figure 9 
provides a homogeneous stable diffusion-controlled flame 2 inches wide 
and up to 30 feet long. Measurements will be made at pressures from .1 
to 2 atm. with temperatures from 1300®K to 3000°K. 

2. CO 2 , CO 

The emissivity of these gases at high temperatures has been 
theoretically calculated by Malkmus [9,10] and by Thomson [11]. For CO 2 
many of these theoretical results have been verified by experiment [ 12 , 
13], These calculations, which are believed to sufficiently represent 
the spectral properties of CO and CO 2 , have been described in detail in 
the original papers and are tabulated in reference 3. For CO, the cal- 
culated values of (S/d)g.pp for the 4. bp. fundamental band are tabulated 
also in reference 3. For CO 2 , the values 


4 

and 


(S/d) 


STP 





(S/d)g,^p 


for the 4.3[a fundamental band and the 2.7|a overtone band have also been 
tabulated, as a function of temperature and wave number. The agreement 
between the theoretical calculations of Malkmus and the experimental 
data of Oppenheim and Ben-Aryeh is reviewed in figure 10. 

3. Carbon Particles 

For LOX-RP engines, the dominant radiating constituent of 
the exhaust plume is carbon particles. For a heat transfer calculation, 
it was necessary to determine the absorption characteristics of this 
particular carbon found in rocket motors. The size of the carbon 
particles must also be known to determine the effect of scattering on 
the heat transfer. Studies were initiated to determine the emission 
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characteristics and size distribution of the carbon particles. The 
spectral emissivity of carbon clouds, produced at the exit plane of 
small rocket motors burning RP-1 and gaseous oxygen, was measured and 
compared with theoretical calculations. The rocket motors for this 
study used Foelsch-type nozzles. These nozzles produced an approxi- 
mately constant temperature and velocity across the exit plane. The 
three motors used were with identical injectors, combustion chambers, 
and nozzle contours as far as the throat section. The diverging sec- 
tions of the nozzles were different and expanded the combustion products 
to area ratios of 5.25, 3.0, and 1.5. In this manner, at a given mixture 
ratio and chamber pressure, carbon formed in a given amount under iden- 
tical conditions was studied at three different temperatures. The pro- 
duct of the carbon mass fraction, Y^, and the carbon cloud’s spectral 
absorption coefficient, c has been determined for wave lengths 
between 1 and 4|a and temperatures between 1000 and 2600°K. The spectral 
dependence of for temperatures below lyOO^'K, is in reasonably good 

agreement with the calculations of Stull and Plass [14] based on the D.C. 
conductivity of bulk graphite and appropriate to very small particles 
(figures 11 and 12). Setting c proportional to and independent 
of temperature is also consistent with the present data below 1750®K. 
However, for temperatures above 17C0°K, c appears to be almost inde- 
pendent of the wavelength. Apparent values of carbon mass fraction, 
determined from the combination of calculation and experiments, decreased 
very strongly with increasing O/F ratios and increased with increasing 
chamber pressure. The carbon particles have also been sampled in these 
small rocket motors and in a full-sized F-1 engine (see figure 18). 

Figure 13 indicates an approximate size of 400 A, which agrees quite 
well with the calculations of Stull and Plass. We thus have the absorp- 
tion coefficients for carbon particle clouds, and since the size of the 
particles is small compared to the wavelengths associated with the 
radiant energy in an exhaust plume, the effect of scattering on the heat 
transfer does not have to be considered. 


VII, RADIANCE CALCULATION FOR INHOMOGENEOUS ROCKET EXHAUSTS 


Two radiative heat transfer computer programs have been written to 
predict radiation from exhaust plumes. These programs differ only in 
geometric considerations, and the less complex of the two, which computes 
radiation from an axisymmetric exhaust plume, is described in this paper. 
The other program is similar but more geometrically complex because it 
was developed to calculate the radiation from a three-dimensional flow 
field such as that encountered with clustered rocket engines. 



The computer programs use a modified Curtis-Godson approximation 
in the well-known radiative heat transfer equation. For purposes of 
calculation, the integral heat transfer equation is replaced by sums, 
over discrete locally homogeneous regions. The band model parameters 
used in these programs are over 25 cm"^ increments . 

The geometry of the axisymmetric program is shown in figure 14. 

The left-hand coordinate system for the plume flow field has the Z-axis 
at the flow field centerline with +Z in the downstream direction. The 
X-Y plane corresponds to the origin of the flow field used (usually the 
nozzle exit), and the X-Y axes are oriented so that the point of interest 
is in the X-Z plane. The point of interest is the origin of a right- 
hand coordinate system with axes U, V, and W. The +W-axis is the outward 
normal to the surface of the unit area receiving radiation, and the U-axis 
is constructed to intersect the Z-axis at some point (this will be oo if 
the axes are parallel). The positive U direction is that given by a 90° 
clockwise rotation of the +W-axis when viewed from Y = -oo. 

Since the U-V plane coincides with the X-Z plane, the location and 
orientation of the point of interest are specified by Z and X coordinates 
and an angle of inclination co. This inclination angle is measured clock- 
wise when viewed from V = -oo between a vector V parallel to and in the 
same direction as the H-Z-axis and the +W-axis. 

The spherical coordinate system used in computing the radiation has 
its origin at the point of interest. The line of sight, or radius vector, 
is and 0 is the angle between i and the 4-W-axis. The angle between 
the projection of ^ in the U-V plane and the 4-U-axis is 0 . 

To compute the radiation flux per unit area, Q/A, the computer pro- 
gram uses constant values of A0, Av, and between the prescribed 
initial and final values in evaluating the summation 


A 



Vj 


0 . V. 

1 1 



o 






) 


sin 0 cos 


A0 A0 Av . 


( 22 ) 


The summation of the transmissivity, t^, begins with = 1 and in 

subsequent increments uses 




exp 



0 



(23) 
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where the optical path, D, taking into account both Doppler- and col 
1 is ion-broadening, is given by 



For collision-broadened lines. 


I 


D 

c 


o 




( 26 ) 


and for Doppler-broadened lines, 


A 


I 


o 




( 27 ) 
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At the center of each ZXg increment along each line of sight, the 
pos ition in the flow field is computed , and the properties at that point 
are determined from the tabulated flow field properties by linear inter- 
polation. This provides an estimate of the pressure, P, temperature, T, 
and concentration, C, at that point, for use in the radiance calculation 
The pressure is converted to a temperature-corrected partial pressure 
(actually a density ratio), using 




p' 

CO 


= -^ (P/1 Am) 


c. 


A table of absorption coefficients listed for the exact values of wave- 
length used in the program and several temperatures is used to determine 
the effective absorption coefficients by linear interpolation in tem- 
perature, The values of 7 , 7 ^, and 1/d are determined using relations 
listed in Table II. Now all or the values are available for evaluating 
a step in the inner summation of equation ( 22 ), and the procedure is 
repeated at the center of succeeding AJl increments until the maximum 
distance, is reached. To save time in searching the flow field, the 

inner sum of equation (22) is simultaneously evaluated for 300 wave 
number values. This is sufficient to cover the spectral range of inter- 
^ est in most problems so that each line of sight has to be traversed only 

one time. 

To account for blockage of radiation along particular lines of 
♦ sight by the engine nozzles and vehicle structure, the program has the 

capability of providing up to 50 "blocking circles." These can be 
located anywhere in the flow field in planes normal to the Z-axis and 
can be either a disk which would block lines of sight passing the plane 
within the circle or a hole which would block lines of sight passing the 
plane outside the circle. Using these blocking circles an engine con- 
tour can be approximated as a stack of "disks," and a structure such as 
an interstage skirt can be approximated as a stack of "holes." 


VIII. APPLICATIONS 


The radiation computer program results are now being compared with 
various radiation measurements; one of these comparisons for which pre- 
liminary results are available is described. The radiation measurement 
was made during a test of a duration model rocket at Cornell Aeronautical 
Laboratory. 
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The short duration rocket testing technique uses nonsteady flow 
techniques with gaseous propellants to provide steady flow model rocket 
firings on the order of 5 to 10 milliseconds. The short test time allows 
high altitude simulation and eliminates the need for nozzle cooling 
systems . 

Short duration models of two F-1 engines (1/45 scale) were used in 
the radiation measurement program. The configuration of the models 
simulated the engine spacing between the center and outboard engines 
on the Saturn S-IC stage. The Warner and Swasey Model 501 fast-scanning 
spectrometer, which was used for making the spectral measurements, is 
capable of making a spectral scan in 1 millisecond with a 0.25 milli- 
second interval between scans. It uses two detectors, and the two chan- 
nels of output are presented on oscilloscopes and photographed. Since 
there is no absolute wavelength reference in the data presentation, 
it must be fixed by the spectral features. Therefore, slight wave- 
length discrepancies in the comparisons presented should be interpreted 
as uncertainty in the data reduction rather than actual comparative 
differences . 

The radiation measurements made in the regions indicated in fig- 
ure 15 have been compared for the F-1 engine at positions 1 and 5. The 
spectrometer field of view at the plume used during the tests was 
approximately 0.8 in. high and 0,25 in. long. No comparisons have been 
made in the impingement region between engines because methods for pre- 
dicting the three-dimensional flow field in this region have not yet 
been perfected. 

The predicted flow field properties and radiation for the F-1 
engine are compared in figures 16 through 19. The comparisons are 
excellent except for a slight wavelength discrepancy in the long wave 
length channel which is attributed to uncertainty in the data reduction 
procedure described previously. 


IX. CONCLUSIONS 


An improved method for calculating radiative heat transfer through 
inhomogeneous, Saturn- type exhaust plumes, has been developed. This 
calculation method uses a band model approach in representing the 

absorption and emission characteristics of the exhaust gases. A modi- « 

fied Curtis-Godson approximation is used in the radiative heat transfer 
calculations which effectively average the band model parameters over 
the inhomogeneous path. This technique has been shown to be accurate 

for strong gradients in temperature and concentration. It has been ^ 

shown that the difference between the measured inhomogeneous trans- 
mittance and that calculated by the modified Curtis-Godson approximation 
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was within the experimental error and was not greater than .02, These 
band model parameters should be completely determined by November 1967. 

An improved set of S/d values have been obtained for H 2 O, CO 2 , and CO. 

A first set of 7 / d values have been obtained for these gases. This 
first set of band model parameters is used in the radiative heat trans- 
fer program. At this time very little experimental data are available 
for comparison with the calculations. As flow field data improve, com- 
parisons between experiment and calculation should be much more favorable 
Extensive flow field work is now in progress. Other improvements will 
come from refined values of the fine structure parameters, particularly 
the dependence of 7 / d on temperature and wave number at long path lengths 
The mass absorption coefficient of carbon as a function of frequency and 
temperature has been determined, using a small perfectly balanced rocket 
nozzle. These data, along with sampling data, show the carbon to agree 
with the theory of Stull and Plass. A three-dimensional method of 
characteristics program and a finite difference technique are being 
developed to predict gas properties in the impinging regions of the 
exhaust plumes of the Saturn V, 






TABLE I I 

RANDOM MODEL 
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R is a normalizing factor 

Xg = S^u/Zwy 
E denotes the mean value 



TABLE III 


Two-Zone Water Vapor Transmittance at 3990 cm 


.1 


Optical Path: 24" per zone 

Ti = T 2 = 1273°K 
= X 2 = 3.49 


Pi (mm) 

PsCmm) 


r- 

Calculated t 



T2 

T 

From 
Eq. (17) 

From 
Eq. (19) 

53 

53 

.720 

.729 

.621 

.628 

.634 

102 

105 

.571 

.566 

.425 

.443 

.450 

150 

153 

.438 

.443 

.284 

.306 

.314 

51 

104 

.743 

.571 

.515 

.524 

.530 

56 

146 

.724 

.464 

.419 

.428 

.435 


TABLE IV 

Two-Zone Water Vapor Transmittance at 3990 cm 
Optical Path: 24" per cone 

Ti = T 2 = 637°K 
Xi = X 2 = 2.53 


Pi (mm) 

P 2 (mm) 

Measured 

Transmittance 

Calculated t 


'^2 

T 

From 
Eq. (17) 

From 
Eq. (19) 

51 

50 

.721 

.728 

.620 

.626 

.634 

102 

100 

.510 

.528 

.368 

.384 

.395 

147 

148 

.374 

.381 

.221 

.241 

.252 

148 

59 

.383 

.677 

.339 

.345 

.355 

105 

49 

.513 

.726 

.462 

.468 

.477 

102 

153 

.507 

.365 

.273 

.286 

.297 


22 

















TABLE V 


Two-Zone Water Vapor Transmittance at 3990 cm"^ 
Optical Path: 8" per zone 
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FIGURE I LOW ALTITUDE CASE — SATURN I EXHAUST PLUME 
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FIGURE 2 HIGH ALTITUDE CASE — SATURN I EXHAUST PLUME 




120,000 FT. ALTITUDE 0/F = 2.25 

OXYGEN/ETHYLENE (CgH^) P. » 1105 PSIA 
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FIGURE 3 !/45 SCALE F-l ENGINES 


DELTA FUNCTION DISTRIBUTION 



distribution functions . 






WAVENUMBER (CM"') 


Figure 6 


Comparison of absorption coefficients ^ obtained by Goldstein 
of the 6.3-p, band at 4T3®K (solid line). Present values are 
given as points . 
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FIGURE 9 END VIEW OF LONG BURNER FLAME 
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FIGURE II. ABSORPTION COEFFICIENT X MASS FRACTION VERSUS 
WAVELENGTH FOR CARBON PARTICLES AT VARIOUS 
TEMPERATURES 
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FIGURE 13 ELECTRON MICROGRAPH COLLECTED SOOT 
PARTICLES-TEST MAGNIFICATION 96500X 
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Properties Predicted for 02/C2H^ Propellants, 0/F = 2.25, 

P( = 1000 psio, Pq‘ 9.61 psfo, and Equilibrium Composition 
During Expansion. 
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Spectral Radiance ( watts / cm^ -/i - ster ) 



Wavelength ( microns ) 

FIG. 17. COMPARISON OF MEASURED AND PREDICTED RADIATION 
NEAR THE EXIT OF A 1/45 SCALE F-1 ENGINE 




Properties Predicted for O 2 /C 2 H 4 Propellents, 
0/F = 2.25, P 5 = 1000 psio, Pj^ = 9.61 psfa , 
and Equilibrium Composition During Exponsion. 
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Remaining Gas Assumed to Be H 2 for Rodiation 
Calculation. No Carbon Wos Predicted Although 
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FIG. 18. PREDICTED PROPERTY VARIATIONS 
9.94 INCHES DOWNSTREAM OF A 1/45 SCALE F-1 ENGINE EXIT 
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